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Abstract 



We analyze the effects of intersite energy correlations on the linear optical 
properties of one-dimensional disordered Frenkel exciton systems. The ab- 
sorption line width and the factor of radiative rate enhancement are studied 
as a function of the correlation length of the disorder. The absorption line 
width monotonously approaches the seeding degree of disorder on increasing 
the correlation length. On the contrary, the factor of radiative rate enhance- 
ment shows a non-monotonous trend, indicating a complicated scenario of the 
exciton localization in correlated systems. The concept of coherently bound 
molecules is exploited to explain the numerical results, showing good agree- 
ment with theory. Some recent experiments are discussed in the light of the 
present theory. 
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I. INTRODUCTION 



Theoretical studies of excited states of quasi one-dimensional (ID) self-assembled molec- 
ular systems, like dipolar coupled molecules in J aggregates,0ll are often based on the well- 
known Frenkel exciton mo del (see Ref. ^ for a comprehensive review). The most remarkable 
optical properties of J aggregates are the narrowing of their absorption spectra (J band)0'i 
and the radiative rate enhancementSi as compared to the monomers, both originating from 
the collectivization of the local molecular states due to the intermolecular dipolar coupling. 
The absorption band is narrowed by about the square root of the so-called number of co- 
herently bound molecules, N*, which, roughly speaking, is the number of molecules covered 
by the exciton wave function.0 The same number determines the factor of radiative rate 
enhancement. A^* is smaller than the physical length of an aggregate since Anderson local- 
ization takes place in ID random systems.! It has been realized that diagonal and off-diagonal 
disorder are particularly important in determining the optical properties of these systems 
because A^* strongly depends on the degree of disorder.0i 

Since the pioneering work by Knapp,i theoretical studies of the optical dynamics of 
Frenkel excitons in ID random systems with intersite correlations have attracted much 
attention.0Hli Until recently, however, it was not possible to experimentally prove the oc- 
currence of intersite correlations. Hence, most theoretical descriptions assumed that disorder 
was indeed uncorrelated. Knoester showed that nonlinear optical techniques are appropriate 
means to determine both the degree of disorder and the correlation length of the diagonal 
disorder.li2l This approach was applied to molecular aggregates of pseudo-isocyanine (PIC) 
by Durrant et al^ and in aggregates of the dye 5,5',6,6'-tetrachloro-l,l'-diethyl-3,3'-di(4- 
sulfobuty)-benzimidazolocarbocianine (TDBC) by Moll et a/.,0 which are known to easily 
form J aggregates. 

One of the general consequences of the occurrence of intersite energy correlations outlined 
in all the studies concerning the exciton absorption band is the increase of its width as 
compared to uncorrelated random systems.B0H3 It was found that the absorption line width 
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approaches the seeding degree of disorder on increasing the correlation length of the intersite 
energy disorder. To the best of our knowledge, however, there was no a detailed study of the 
factor of radiative rate enhancement, carrying out information about the coherently bound 
molecules. In a recent paper devoted to pairwise correlated diagonal disorder (site energies 
appear at random to pairs of nearest-neighbor sites) ,0 we found the counterintuitive result 
that this factor drops upon the occurrence of correlations. This means that excitons at 
the bottom of the band become more localized after introducing pairwise correlations while 
the underlying lattice is less disordered, in the sense that only half of the site energies are 
truly random variables. Stronger localization near the band edge upon introducing the 
correlations in the disorder was also pointed out in the recent paper by Russ et a/.0 

In this paper, a detailed study of ID molecular systems with intersite energy correla- 
tions is performed in order to uncover the effects of the short-range to long-range correlation 
length crossover on the linear optical response of Frenkel excitons. We show that the factor 
of radiative rate enhancement drops upon increasing the energy correlation length whenever 
this length is smaller than a given characteristic length. For larger energy correlation lengths 
the behavior of the factor of radiative rate enhancement is the opposite and increases on 
increasing the correlation length. We provide a simple recipe for estimating this character- 
istic length. The absorption band width presents a monotonous increase upon rising the 
energy correlation length and approaches the width of the distribution of seeding disorder 
when the correlation length exceeds the above mentioned characteristic length. We also use 
the concept of coherently bound molecules to interpret numerical data. Applying this con- 
cept in a self-consistent way, we find the scaling of the factor of radiative rate enhancement 
and the width of exciton absorption line with the energy correlation length and show that 
theoretical estimates reasonably fit the results of numerical simulations. The remainder of 
the paper is organized as follows. In the next Section we formulate the model we will be 
dealing with. Section |IT1| deals with the theoretical estimates which are then used in Sec. ^ 
for explaining the resutls of numerical simulations. Some recent experiments are discussed 
in the light of the present theory in Sec. 0. Section ^ summarizes the paper. 



II. MODEL OF DISORDERED SYSTEM 



We consider N :$> 1 optically active, two-level molecules forming a regular ID lattice with 
spacing unity. To build up our correlated disordered lattice with intersite energy correlation 
length Nc, we create S = N/Nc consecutive segments with equal transition energies 
within each segment. The set {es}f=i ^-re statistically independent Gaussian variables with 
probability distribution 




and (e^es') = cf'^Sss', where brackets denote the average over the join energy distribution 

p{{^s]) = i{p{^s). (2) 

and a is the seeding degree of disorder. 

For our present purposes, we will neglect all thermal degrees of freedom (electron-phonon 
coupling and local lattice distortions) and assume that disorder originates from a static offset 
of the on-site energies (diagonal disorder). The effective tight- binding Hamiltonian for the 
ID Frenkel-exciton problem in the random lattice with correlated disorder can be then 
written as follows 

S Nc N 

n^J2^sJ2\^ + Ncis-l)){n + Nc{s-l)\+ Jnm\n){m\, (3) 

s=l n=l n,m=l 

where {u) is the state vector of the i^th molecule and Jnm — J\n-m\ are the hopping integrals, 
which are assumed to be negative, as it takes place for J aggregates, and not subjected to 
disorder. 

We are interested in the absorption line shape and its width as well as in the factor of 
enhancement of the exciton radiative rate relative to the monomer spontaneous emission 
rate. The absorption line shape is calculated as 

1 / ^ \ 
m = j;^{i:^^jS{E-E,)y (4) 
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Here we assume that the system size is much smaller than the optical wavelength. The 
index j runs over all eigenstates of the Hamiltonian (|^) and Ej denotes their corresponding 
eigenenergies. The oscillator strength of the j th eigenstate with components a^"^ is given 
by 

- (E -f) , (5) 

where the dipole moment of an isolated monomer is taken to be unity. 

The factor of radiative rate enhancement is defined through the average oscillator 
strength per state at energy idUli 

AiE) ^ (6) 
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where 

PiE) = j;^{T.HE~E,)^ (7) 

is the normalized density of states. We take max{/i^y(E)} as a measure for the enhance- 
ment of the exciton radiative ratei or, in other words, for the number of coherently bound 
molecules A^*.i In what follows, the above magnitudes will be determined both analytically 
and numerically as a function of the correlation length Nc- 



III. THEORETICAL ESTIMATES 

Having presented our model, we now describe the method we have used to estimate the 
factor of radiative rate enhancement A^* and the linear absorption band width ai. To this 
end, we rewrite the Hamiltonian (^) in the excitonic representation using the eigenfunctions 
of the unperturbed Hamiltonian (with = for s = 1, . . . , S). For the sake of simplicity, 
we assume here periodic boundary conditions. The Bloch plane waves are then the exact 
eigenfunctions of (j^) in the absence of disorder: 



where K = 2Tik/N belongs to the first Brillouin zone (0 < < A^). In the /^-representation, 
the Hamiltonian (^) reads 

n = Y.EK\K){K\ + ^KK'\K){K'\, (9a) 

K K,K' 

where Ek is the unperturbed exciton energy 

N 

Ek = T. (9b) 

n=l 

and Akk' is the intermode scattering matrix 

1 S Nc 

^KK' = ^ E E e^(^-^')['^+^=(^-i)l. (9c) 

s=l n=l 

The role of Akk' is two-fold. The off-diagonal elements Akk' {K ^ K') mix the exciton 
states. If the typical fluctuation of Akk' does not exceed the energy difference between 
the states K and K' (the so called perturbative limit) then the energy shifts of the corre- 
sponding states, given by the diagonal elements Akk^ is the main effect of the disorder and 
result in inhomogeneous broadening of exciton levels for an ensemble of chains. The typical 
fluctuation of A^k has a direct relationship with the inhomogeneous width of the exciton 
state \K). 

For nonperturbative magnitudes of the disorder, off-diagonal elements of the scattering 
matrix, A^k' {K 7^ K'), mix the exciton states, resulting in their localization on chain 
segments of a typical size smaller than the chain length and subsequently affecting the 
exciton optical response. Recall that, for a perfect circular chain, only the state \K = 0) 
is coupled to the light and carries the entire exciton oscillator strength, which is then 
times larger than that for an isolated molecule. Being mixed with other (nonradiative) 
states {K 7^ 0), the radiative state loses a part of the oscillator strength due to its spreading 
over the nonradiative ones. Thus, the number of coherently bound molecules A^* < 
arises as the enhancement factor of the oscillator strength of the localized exciton states.0 
Accordingly, the inhomogeneous width of the optical exciton line will also be subjected to 
renormalization since A^ should be replaced by A^*.i 
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A. Perturbative treatment 



Assuming the perturbative limit (a quantitative condition of its validity will be given 
below), let us gain insight into the magnitude of the typical fluctuation of the scattering 
matrix IS.kk'- To this end, one should calculate its second moment using the joint distribu- 
tion (^) of the seeding fluctuations (being expressed trough linear combinations of Gaussian 
variables with zero mean, the Akk' are also Gaussian variables with {Akk') = 0). The 
magnitude of our interest will be cr^x' = (l^-ft'X'P)- It is an easy task to arrive at 

sin^ [{K - K')N J 2] 



(Jr. 



(10) 



NN, sm''[{K - K')/2] ' 

In the limits K, K' <^ 1 and {K — K')Nc -C 1, which are of major importance from the 
viewpoint of optical properties, Eq. (|T0|) reduces to 



This expression generalizes the result we obtained previously for pairwise correlated disorder 
{Nc = 2)lll to the case of an arbitrary energy correlation length N^. As follows from Eq. (pT]), 
the motiona narrowing effect (A^~^ scaling) is weakened here by the factor and is deter- 
mined now by the chain length counted in units of the energy correlation length Nc or, 
in other words, by the number of correlation segments 5* = N/N^. The motional narrowing 
effect completely disappears at = N. Similar conclusions were drawn in Refs. [7|,PHn]. 

As it was already mentioned above, the state \K = 0) carries the entire oscillator strength 
of the system so that the optical absorption spectrum presents an isolated Gaussian peak 
centered around the eigenenergy Ek=o and it is characterized by its standard deviation 

-1 = • (12) 

This result holds whenever the state K = does not become mixed with the nearest one \Ki), 
where Ki = 2tt/N. The mixing is governed by the off-diagonal element A^io- Hence, one 
should compare axio with the energy difference Ek^ — Eq. The perturbative approach is still 
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applicable provided axio < Ek^ —Eq and fails otherwise. Thus, the equality (TkiO = Ek^—Eq 
determines a value of a (for fixed both chain length and energy correlation length N^) 
which separates the ranges of perturbative and nonperturbative magnitudes of disorder. 
Taking Ek in the nearest-neighbor approximation 

= -2 J cos (fsT) ^ -2 J + JJsT^ (13) 

where J is the nearest-neighbor coupling and is assumed to be large, one gets 

^ 47rV 

which determines the crossover between the perturbative and non-perturbative limits. 

B. Nonperturbative limit 

At nonperturbative magnitudes of disorder, not all molecules of the chain contribute to 
the optical spectra of the whole ensemble, but only a portion of them being a function of 
the degree of disorder cr/J.i The reason is the localization of the excitonic states arising 
from disorder. Therefore, the number of coherently bound molecules, A^*, (averaged number 
of molecules covered by optically active localized exciton states) should replace the number 
of molecules in the system, A^. In Ref. |19|, a simple rule for estimating A^* was formulated 
whenever the disorder is uncorrelated. It exploits the findings that i) the lowest localized 
exciton states can be classified in several groups of states (two or sometimes three), each 
one localized on a certain chain segment of a typical size A^*, ii) each segment does not 
overlap with the other^ll (see also Refs. pO| , pT| and Sec. |^ for more details) and, what i 
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most important, the states of each group have the energy structure similar to that for an 
homogeneous chain of size A^*, i.e., approximately given by Eq. (|13|) with A^ replaced by A^*. 

The rule proposed in Ref. |19| for estimating A^* consists simply of applying the pertur- 
bative criterion (|1^) to a typical localization segment, i.e., substituting A^ by A^* in Eq. ( p!^ 



and considering now A^* as an unknown parameter. As was shown in Refs. 0,^, this ap- 
proach works surprisingly well in fitting the numerical data concerning the optical response 



of ID Frenkel excitons. In doing so, one arrives at the following estimate of the number of 
coherently bound molecules 

According to Ref. the standard deviation of the exciton absorption spectrum cti can be 
estimated by Eq. (|l^), again replacing by A^*, that gives us 

The Ac scaling laws of A^* and ai are the main results of the present paper. Note that 
cr-scaling of A^* and o"i was previously reported in Refs. ^|19] and |9| JT8| , p!9|j2^ -p^, respectively. 



We stress that the number of coherently bound molecules drops while the standard deviation 
of the exciton absorption spectrum rises upon increasing the intersite energy correlations, 
meaning stronger scattering of the exciton in a more correlated system. This somewhat 
unexpected result stems from the fact that motion narrowing effect is partially suppressed 
as intersite energy correlations occurs, yielding a smaller reduction of the seeding disorder 
as compared to the case of uncorrelated disorder. 



C. Short- to long-range correlations crossover 

The concepts presented in the previous subsection as well as the scaling laws (|15D and (|16D 
are valid within the range Nc < N*. However, both numbers, Nc and A^*, will approach 
each other as the correlation length Nc rises and, starting from a certain value of Nc we 
will obtain A^* ^ Nc- We now provide arguments how to estimate the magnitude of Nc at 
such a crossover. As soon as A^* and Ac approach each other, there is no reason for the 
reduction of disorder due to the motional narrowing since the energy of any segment of size 
Ac is fixed, i.e, does not fluctuate. Thus, the typical energy mismatch between two adjacent 
segments is of the order of a. On the other hand, the fact that A^* ^ A"c means that the 
coupling between adjacent segments is of the order or smaller than their energy separation 
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(T. From this it follows that the condition we are looking for is precisely that under which 
two adjacent segments of size A^^^ start to be decoupled of each other. 

In line with this reasoning, let us consider the two segment problem, taking for the sake 
of simplicity the corresponding Hamiltonian in the nearest-neighbor approximation 

r Nc~i 



ea\<y,n){n,a\ — J ^ (\a,n){n + l,a\ + \a,n + l){n, 



a=a,b Ln=l n=l 



- j(|a,iVe)(l,6| + |fe,l)(iV„a|), (17) 

where a and h are the segment labels and the second term describes the coupling of the 
segments. Making use of the excitonic transformation 



K 



where K = nk/ {N^ + 1) and k ranges from 1 to N^, one passes from (0) to 

n= E T.i^a + EK)\a,K){K,a\+ Y.i^ah)KK'\a,K){K' ,h\ . (19a) 

a=a,b K K,K' 

Here, Ek is given by Eq. (|T3D and the interaction Hamiltonian reads 

9 T 

{nab)KK' = (-i)'^']^rTT ""^^^^ ""^^^'^ ■ ^^^^^ 

The interaction Hamiltonian {Ti.ah)KK' couples excitonic modes belonging to different seg- 
ments and thus leads to their spreading over both segments. Recall that, from the viewpoint 
of the optical response, mixing of the lowest modes with Ki = K[ = -k /{N^. + 1) is of major 
importance. We aim to elucidate for which segment size Nc this coupling can be considered 
as a perturbation. Remembering that the typical energy mismatch between the segments 
is ~ cr and taking \{'Hab)K^K[\ = c", one arrives at the following estimate of the size we are 



looking for 



iv: + i = f^l , (20) 



where the condition Ai'c ^ 1 was used. 
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Summarizing, the typical energy separation of two adjacent segments starts to dominate 
over their couphng at Nc > N* so that the segments can be considered as decoupled of each 
other. As a consequence, the localization length of the lowest exciton mode will be the same 
that the energy correlation length Nc- 



IV. DISCUSSION OF THE NUMERICAL RESULTS 

In order to calculate the quantities of interest defined by Eqs. (0)-(0) we replaced 6{E — 
Ej) by the Heaviside step-function {1/R)6{\E — Ej\ — R), where R is the resolution, taken to 
be 0.05 in all the calculations. For energy units, J = 1 was chosen. The standard deviation 
a was also fixed to be 0.2. The only varied parameter was the energy correlation length Nc- 

We diagonalized the Hamiltonian (Q) for open chains of size = SN^, where for a given 
Ac; the number 5* was chosen such that A^ would be the closest one to 250. The number 
of randomly generated systems was 1000 for each value of A^. We do not present here the 
exciton absorption spectra since they are rather standard ones, showing the characteristic 
asymmetry and red shift due to disorder discussed in detail for uncorrelated disorder in 
Refs. |9| JT8|j25| and for the correlated case in Refs. [7|, [TllJl3| . In Fig. ||, the standard deviation 



of the absorption spectra ai is shown as a function of the energy correlation length A^. To 
accurately determine the standard deviation ai we fitted the low-energy side of the spectra 
using Gaussians. One can see from Fig. |l| that ai monotonously increases as Nc rises 
and finally approaches the value 0.2 (= a) for the case at hand. The crossover, occurring 
approximately at A^* = 5, is in good agreement with the estimate A"* ^ 4 given by Eq. (^0|). 

Figure |^ shows the behavior of the factor of radiative rate enhancement max{yU^y(E)} 
as a function of A^c- First, it drops upon increasing Nc and then it grows for A^c > 4, also 
in full correspondence with our theoretical prediction done in Sec. |ITT[ Strong fluctuations 
observed in Fig. ^ result from the fluctuations in the density of states: the latter appears in 
the denominator in the definition of max{/i^^(E)} [see Eq. (^)]. We do not present in Fig. |^ 
how the quantity max{/x^y(E)} behaves for A'c > 15 for several reasons. First, as soon as the 
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segments of correlated energies start to be decoupled of each other (the case for Nc > 4, see 

Fig. ^ and the oscillator strength per state fJ'l^iE) behaves monotonously, i.e., has no 

maximum (see Appendix thus proving to be an inconvenient measure of the extension 

of the optically active exciton states. Second, it is evident that, for decoupled segments, the 

localization length A^* is equal to the energy correlation length Nc- Figure |^ confirms the 

last statement. There, we presented typical realizations of the eight lowest eigenfunctions 

of the Hamiltonian (|^) for three values of the energy correlation length Nc = 4, 10, and 40. 

They demonstrate that the energy segments actually turn out to be decoupled of each other 

as the segment size rises. Indeed, for Nc = 4, the typical extension of the eigenfunctions, 

A^*, of the lowest excitonic states is more than 10 (the oscillator strengths /i^ of nodeless 

states carry this information). Note that it sufficiently exceeds the energy correlation length. 

For larger value, A^ = 10, both magnitudes approach each other, although still keeping the 

relation A^* > Ac. Besides, in the last case, the eigenfunctions tend to be grouped on 

different segments of size ~ A^. Finally, for the largest presented value of Nc{= 40), the 

eigenfunctions fKn, oscillator strengths and eigenenergies within each group follow 

fairly well the corresponding formulae for an isolated segment of size A^: 

/ 2 

'^Kn = j,j , , sin (Kn) , (21a) 



Nc + l 



1 



1 - (-1)' 



^cot^(—), (21b) 



2{Nc + 1 

SK = e + Ek, (21c) 

where e is the site energy of the molecules within each segment. Note that such a structure 
(called hidden in Ref. |1^) of the exciton low-energy spectrum is also present for A"c = 4 but 
scaled in average by the number A^* instead of Nc- 

V. DISCUSSION OF THE EXPERIMENTAL RESULTS 

In the recent experimental works by Durrant et al^ and by Moll et al^ it was reported 
on the presence of strong correlations in diagonal disorder in molecular aggregates of PIC 
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and TDBC, respectively. Moreover, it was found in both studies that the energy correlation 
length Nc is of the order of the localization lengthlll and even larger.lll This limit strongly 
suggests that almost no motion narrowing effect has to be detected.0'0 Nevertheless, the 
observed J band width in both cases was definitely motionally narrowed. The possible 
ordering of the first solvatation shell in the presence of highly ordered aggregate during the 
freezing process and the subsequent reduction of the local disorder was given as a possible 
explanation for the J band narrowing. 

It seems to be interesting to estimate the critical length A*"* for the above aggregates 
on the basis of the experimental data presented in Refs. |15|,|16[ Taking J = 600 cm~^ and 
a = 13.5 cm^^ for J aggregates of PIC,0 one gets A*"* ^ 9 what is consistent with the 
relation Nc — N* found by Durrant et al^ For J aggregates of TDBC (J = 850 cm^^ and 
a = 67 cm^^), one arrives at A^* ~ 6. Thus, the experimental estimation of the energy 



correlation length Nc to be of the order of several hundreds molecules obtained in Ref. |6 
contradicts, in our opinion, the magnitude of localization length N* ^ 40 extracted from 
the spectroscopic data. 



VI. SUMMARY AND CONCLUDING REMARKS 

In summary, we have studied the effects of intersite energy correlations on the character- 
istics of linear absorption spectrum corresponding to a weakly localized ID Frenkel exciton. 
We have found that the factor of radiative rate enhancement shows a nonmonotonous behav- 
ior upon rising first decreasing as A^^"^^^ ^^'^ then growing linearly. The crossover occurs 
when the typical coupling of segments is smaller than the typical separation between their 
energy levels, being of the order of the width of the seeding Gaussian energy distribution a. 
The linear trend indicates that the segments become independent one from the other when 
Nc > (27r^J/o")^/^. The standard deviation of the absorption line monotonously increases as 
Np for low energy correlation length and saturates when Nc > (27r^ J/a)^/'^, namely when 
the localization length A^* coincides with the energy correlation length A^^. Our results lead 
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us to a better understanding of recent experimental results in PIC and TDBC molecular 
aggregates. 
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APPENDIX A: EVALUATION OF fil,iE) FOR DECOUPLED SEGMENTS 



Let us calculate nJ^^E) = p{E)/I{E), assuming the energy segments to be decoupled of 
each other. Under such an assumption, one can rewrite p{E) and I{E) as follows 



S Nc 



EE HE-es-EK) 



s=l k=l 



(Ala) 



1 / 



S Nc 



m = J^{J:Y.^^mE-es-EK) 



(Alb) 



s=l k=l 

where Ek = — 2Jcos[7rA;/(A^ + 1)], k = 1,2, Nc and {e^jf^i is a stochastic Gaussian 
energy offset distributed according to Eq. (^. Making use of the integral representation of 
the (5-function 6{x) = (l/27r) dt exp{ixt) , it is easy to carry out the average in Eqs. ([Ala]) 
and ( |Alb| ) explicitly to arrive at 



1 ^'^ / 1 



1/2 



exp 



{E - Ek) 
2ct2 



(A2a) 



Nc 



m = - E /^^ 



k=l 



1/2 

2 , exp 



{E-E_ 



2a2 



(A2b) 



Accounting further for the fact that the lowest exciton state (with Kx = ii/[Nc + 1]) carries 
almost the entire oscillator strength J2k f^K = 8/^^; namely 8%) and thus keeping in 



Eq. ( A2b ) only the first term, we finally get 

^2 r ^E^-E 



y^-J{E) 



TT 



IN,., 



1 + E exp 

k=2 



Ki { _ Ek + Ek^ 



(A3) 
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From this it follows that fi~^{E) monotonously rises from the value n'^/SNc to infinity, when 
passing from £' = —00 to £' = 00. As a consequence, fil^{E) monotonously decreases from 
the value SNc/tt"^ to zero. It should be pointed out that the characteristic energy scale of 
the exponentials in Eq. ( [A3| ) is determined by the ratio {Ek — EKi)/cr'^ which goes down as 
the segment size Nc rises. This means that close to the band bottom, which is precisely the 
energy region we are interested in, the contribution of the exponentials in Eq. ( |A3|) may be 
large. Then, the magnitudes of H%j{E) obtained from numerical simulations will strongly 
differ from the asymptotic value SNc/tt"^. Due to the peculiarity outlined above, the oscillator 
strength per state fil^{E) appears to be an inconvenient measure for the extension of the 
optically active exciton states for the case of larger energy correlation length, when the 
energy segments are decoupled one from the other. 
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FIGURES 

FIG. 1. Standard deviation of the exciton absorption spectra ai versus the energy correlation 

length Nc- The degree of seeding disorder is c = 0.2. 

FIG. 2. Factor of the radiative rate enhancement, max{//^y(E)}, versus the energy correlation 
length Nc- The degree of seeding disorder is a = 0.2. 

FIG. 3. Typical realizations of the lower exciton eigenmodes for three values of the energy 
correlation length a) Nc = 4, b) 10, and c) 40. 
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